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Abstract 

We compute the branching fractions of B s d — * in the type-II two-Higgs-doublet 
model with large tan (3. We find that the parameters of the neutral Higgs sector of the 
two-Higgs-doublet model cancel in the result, so that the branching fractions depend 
only on the charged Higgs mass and tan (3. For large values of tan (3 and a charged 
Higgs mass above the bound from b — > sj, we find that the branching fractions can 
be enhanced by up to an order of magnitude or suppressed by up to a factor of two 
compared to the Standard Model result. We point out that previous calculations in the 
literature are gauge-dependent due to the omission of an important diagram, which gives 
the dominant contribution in the 't Hooft-Feynman gauge. We have analyzed in detail 
the region of the (M H +, tan (3) plane to be probed by searches for B s — > in Run 

II of the Tevatron. Since the branching fraction increases like tan 4 (3, this decay mode is 
complementary to b — > 57 and efficiently probes the large tan (3 region. For tan (3 = 60, 
an integrated luminosity of 20 fb -1 in an extended Run II will probe charged Higgs 
masses up to 260 GeV, if the background to B s — > is small. For the same value 

of tan (3, the LHC may be able to explore charged Higgs masses up to 1 TeV using this 
decay. 

PACS: 13.20.He, 12.60.Fr, 14.80.Cp 
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1 Introduction 

The ongoing and forthcoming high-statistics B-physics experiments at BaBar, BELLE, 
HERA-B, the Tevatron, and the LHC experiments ATLAS, CMS and LHCb [Jl|] will 
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probe the flavor sector of the Standard Model (SM) with high precision. These experi- 
ments may reveal physics beyond the SM, and a substantial theoretical effort is devoted 
to calculating the observables that will be tested in various scenarios of new physics. 

A common feature of all popular weakly-coupled extensions of the SM is an enlarged 
Higgs sector. In this paper we study the type-II two-Higgs-doublet model (2HDM), 
which has the same particle content and tree-level Yukawa couplings as the Higgs sector 
of the Minimal Supersymmetric Model (MSSM). If the ratio tan/? of the two Higgs 
vacuum expectation values is large, the Yukawa coupling to b quarks is of order one and 
large effects on B decays are possible. Direct searches for the lightest neutral MSSM 
Higgs particle have begun to constrain the low tan (3 region in the MSSM, because the 
theoretically predicted mass range increases with tan (3. Hence observables allowing us 
to study the complementary region of large tan (3 are increasingly interesting. A further 
theoretical motivation to study the large tan (3 case is SO(10) grand unified theories [§]: 
they unify the top and bottom Yukawa couplings at high energies, corresponding to 
tan (3 of order 50. 

The leptonic decay B$ — > £ + £~, where d! — d or s and £ = e, fi or r, is especially 
well suited to the study of an enlarged Higgs sector with large tan (3. In the SM the 
decay amplitude suffers from a helicity-suppression factor of mi/nib, which is absent in 
the Higgs-mediated contribution. This helicity suppression factor numerically competes 
with the suppression factor of {mi/M w ) tan (3 stemming from the Higgs Yukawa cou- 
plings to the final state leptons, so that one expects the new contributions in the 2HDM 
to be similar in size to those of the SM. 

Earlier papers have already addressed the decay B d i — > £ + £~ in the 2HDM or the full 
MSSM Yet the presented results differ analytically and numerically substantially 

from each other, so that we have decided to perform a new analysis. 

This paper is organized as follows. In section ^| we give a brief review of the type-II 
2HDM. In section |3] we review the SM calculation of the decay Bd> — > £ + £~ and de- 
scribe our calculation of the relevant 2HDM diagrams. We finish section |3] by combining 
the results for the 2HDM diagrams and giving compact expressions for the branching 
fractions. In section |] we compare our result with previous calculations. In section^ we 
present a numerical analysis of our result and estimate the reach of future experiments. 
We present our conclusions in section Finally the appendix contains a discussion of 
trilinear Higgs couplings. 



2 The two-Higgs-doublet model 



3 



2 The two-Higgs-doublet model 

In this paper we study the type-II 2HDM. The model is reviewed in detail in ref. P]. 
The 2HDM contains two complex SU(2) doublet scalar fields, 

*-(«)■ MS)- 

which acquire the vacuum expectation values (vevs) (0°) = Vi and break the elec- 
troweak symmetry. The Higgs vevs V\ and t> 2 are constrained by the W boson mass, 
M^r = ~# 2 (ui + uf) = \g 2 v 2 SM , where u SM = 174 GeV is the SM Higgs vev. Their 
ratio is parameterized by tan (3 = v 2/v 1. 

Since in this paper we are not interested in CP-violating quantities, we assume CP is 
conserved by the Higgs sector for simplicity. The mass eigenstates are then given as 
follows. The charged Higgs states are 

G + = <f>t cos (3 + <\>X sin p 

H + = -<j>ffanP + <pZcaaP, (2) 
and their hermitian conjugates. The CP-odd states are 

G° = <pi' 1 cos (3 + 0°'* sm P 

A = -^/sinp + ^cosp, (3) 

where we use the notation 0° = v i + ^(0i' r + ^0°'*) for the real and imaginary parts of 
0°. The would-be Goldstone bosons and G° are eaten by the W and Z bosons. The 
CP-even states mix by an angle a giving two states, 

uO /0,r . ,0,r ■ 

H = (pi cos a + 02 sm a 

h° = — (j)i' r sma + 02' r cosa. (4) 

In order to avoid large flavor-changing neutral Higgs interactions we require natural 
flavor conservation [JT0|] . We impose the discrete symmetry $1 — > — $1, $ 2 — >■ $2 
(which is softly broken by dimension-two terms in the Higgs potential), with the SU(2) 
singlet fermion fields transforming as d — > —d, u — > u, e — > — e. These transformation 
rules define the type-II 2HDM and determine the Higgs-fermion Yukawa couplings. The 
Yukawa terms in the Lagrangian are: 



£ Yuk = -Y d Q^d - Y u Q$ c 2 u - YiL^e + h.c. 



(5) 
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where $ c = ir 2 &*. Down-type quarks and charged leptons (up-type quarks) are given 
mass by their couplings to $ x ($2). 

Most of the Higgs couplings needed in our calculation are given in ref. [j^]. In addition 
we must consider the trilinear H + H H couplings (H = h°, H°) which were first given 
in ref. [^J and are discussed in appendix 



3 Effective hamiltonian for B — > £ + £ 

The decay B d / — > £ + £~ proceeds through loop diagrams and is of fourth order in the 
weak coupling. In both the SM and 2HDM, the contributions with a top quark in the 
loop are dominant, so that one may describe the decay at low energies of order m b by a 
local bd'££ coupling via the effective hamiltonian, 

H = % aEA < & [CsQs + C P Qp + C A Q A ] . (6) 
V2 27r sin W 

Here Gf is the Fermi constant, a em is the electromagnetic fine structure constant and 
Ow is the Weinberg angle. The CKM elements are contained in £ t = V t * b V t d>- The 
operators in (g) are 

Qs = m b bP L d'££ } Q P = m b bP L d'l 75 £ } Q A = bfP L d! hffct, (V 

where Pi = (1 — 75) /2 is the left-handed projection operator. We have neglected the 
right-handed scalar quark operators because they give contributions only proportional 
to the d! mass. The vector leptonic operator Iry^i does not contribute to — > £ + £~, 
because it gives zero when contracted with the B# momentum. Finally, no operators 
involving a^ v — i [7^, j„]/2 contribute to B d t — > £ + £~. 

Because mj <C Mw,fn t , M H +, there are highly separated mass scales in the decay 
Bd' — ► £ + £~ ■ Short-distance QCD corrections can therefore contain large logarithms 
like \og(m b /M w ), which must be summed to all orders in perturbation theory with 
the help of renormalization group techniques. The calculation of the diagrams in the 
full high-energy theory gives the initial condition for the Wilson coefficients at a high 
renormalization scale u on the order of the heavy masses in the loops. The hadronic 
matrix elements, however, are calculated at a low scale [i = 0(rrib) characteristic of 
the B d / decay. The renormalization group evolution down to this low scale requires the 
solution of the renormalization group equations of the operators Qa, Q s an d Q p- Yet the 
operator Qa has zero anomalous dimension because it is a (V — A) quark current, which 
is conserved in the limit of vanishing quark masses. Similarly, the operators Qs and Q p 
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Figure 1: Dominant SM diagrams. 

have zero anomalous dimension because the anomalous dimensions of the quark mass 
m&(//) and the (chiral) scalar current bP^d'^) add to zero. Hence the renormalization 
group evolution is trivial: if the bottom quark mass in Q s and Qp is normalized at a low 
scale \i = 0(rrib), then no large logarithms appear in the effective hamiltonian or in the 
decay rate. 

In the SM, the dominant contributions to this decay come from the W box and Z pen- 
guin diagrams shown in fig. [I]. These diagrams were first calculated in JTT| ] and give a 
non-negligible contribution only to the Wilson coefficient Ca- There is no contribution 
from a photonic penguin because of the photon's pure vector coupling to leptons. There 
are also contributions to the Wilson coefficient C$ from a SM Higgs penguin [ |12~| ] and 
to the Wilson coefficient Cp from the would-be neutral Goldstone boson penguin [JT3|] , 
but these contributions to the amplitude are suppressed by a factor of m^jM^ relative 
to the dominant contributions and can be ignored. 
The SM decay amplitude is then given by the Wilson coefficient [ |TT| ] 

C A = 2Y(x t ), (8) 

where x t = m^(mt)/M^ = 4.27 ± 0.26 and m t is evaluated in the MS scheme at 
H = m t , giving m t (m t ) = 166 GeV. The function Y(x t ) is given by Y(x t ) = Y (x t ) + 



6 



B Sy d in a Two-Higgs-doublet Model 



6 d' 



W 



H+ 



h°H°,A° 



h°H°,A° 



I I 
b d' 



H+,G + 



h°,H°,A° 



Figure 2: Dominant diagrams in the 2HDM with large tan (3. 



2*-Yi(x t ), where Y (x t ) gives the leading order (LO) contribution calculated in [ JTTj ] and 
Yi(x t ) incorporates the next-to-leading (NLO) QCD corrections and is given in [JT3]]. 
The NLO corrections increase Y(x t ) by about 3%, if m t is normalized at [i = m t . 
Numerically, 

1.55 



Y(x t ) = 0.997 



m t [m t ) 
166 GeV 



(9) 



where we have parameterized the dependence on the running top quark mass in the MS 
scheme. 

We limit our consideration to the case of large tan (3, for which the 2HDM contributions 
to this decay are significant. In the large tan/3 limit, the Wilson coefficients Cp and 
Cs receive sizeable contributions from the box diagram involving W and H + and the 
penguins and fermion self-energy diagrams with neutral Higgs boson exchange shown 
in fig. ^ There are no new contributions to Ca in the 2HDM, which therefore retains its 
SM value. 

We have calculated the individual diagrams explicitly in a general gauge, keeping 
only the terms proportional to tan 2 [3. Although each diagram that involves a W ± or 
G ± boson is gauge-dependent, their sum is gauge-independent. For compactness, we 
give the results of the individual diagrams below in the 't Hooft-Feynman gauge. 
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3.1 Box diagram 

The box diagram in fig. Ogives the following contribution to Cs and Cp\ 

box m e .„„2 



C 



box 



c 



2M 2 ? 



tan p B + (x H +, x t ) 



(10) 



Here xh+ = Mjj+ /M^ and x t was defined after (]8|) in terms of mt(m t ). Strictly speak- 
ing, in a LO calculation like ours, one is not sensitive to the renormalization scheme 
and we could equally well use the top quark pole mass in x t . However, experience with 
NLO calculations in the SM [[T3]] shows that with the definition of m t adopted here, 
higher-order QCD corrections are small in leptonic decays. Finally the loop function 
B + in O reads 

log y log x 



B+(x,y) 



y 

x-y 



y-i 



X 



(11) 



B + (xh+, x t ) also contains the contribution from internal up and charm quarks with 
m c = m u = from the implementation of the GIM mechanism. The effect of a nonzero 
charm quark mass is negligibly small. 



3.2 Penguins 

The penguin diagram with H + and W + in the loop (see fig. □) contributes 



u 5 



— tan f3P + (x H+1 x t ) 
tan 2 (3 P + (x H +,x t ) 



sin 2 a 



+ 



cos 2 a 
M 2 HQ 



Ml 



(12) 



Here again all three quark flavors enter the result from the GIM mechanism, and the 
effect of nonzero charm quark mass is negligible. By contrast, in the penguin diagram 
involving H + and G + in the loop only the internal top quark contribution is relevant, 
because the coupling of G + to quarks is proportional to either of the quark masses and 
we neglect m s . This diagram gives 



n peng,2 



^jpeng,2 



17l£ 2 

— tan pP + (x H +,x t ) 



tan 2 f3 P + (x H +,x t ) 



sin 2 a (M 2 H+ - Ml) 



Ml 



+ 



cos 2 a(M 2 H+ -M 2 m/ 



M 2 H0 
M 2 H+ 



M 2 o 



M 2 W 
M% 



M 2 W 



(13) 
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The results in (jT^) and ( Jill ) involve the loop function 



a; - ?/ 



x log a; y logy 
.t — 1 y — 1 



(14) 



3.3 Self-energies 

Before we write down the result from the diagrams with self-energies in the external 
quark lines, we discuss how these contributions come into play. A treatment of flavor- 
changing self-energies has been proposed in [ |l5| , |T6] ] and analyzed in some detail in JT7[| . 



In these works flavor-changing self-energies have been discussed in a different context, 



the renormalization of the H^-boson coupling to quarks. In [|I5|-|T7|] counterterms have 
been chosen in such a way that the flavor-changing self-energies vanish if either of the 
involved quarks is on-shell. For the down-type quarks these counterterms form two 
3x3 matrices in flavor space, one for the left-handed quark fields and one for the right- 
handed ones, and similarly for the up-type quarks. It was argued in [ |16] , |T7l ] that the anti- 
hermitian parts of the counterterm matrices for the left-handed fields can be absorbed 
into a renormalization of the CKM matrix, and the hermitian parts of the matrices can 
be interpreted as wave function renormalization matrices Z^ and Zf- with i,j = d,s,b 
for our case of external down-type quarks. 

However, it has also been argued Jl8| ] that the on-shell scheme of [|i6jp7l| is not gauge 
invariant. In addition we find that the approach of [0-f7J| leads to an inconsistency 
in our calculation. We cannot cancel the anti-hermitian parts of the self-energies in the 
external lines with the counterterms for the CKM matrix, because unlike in the case 
of the W coupling renormalization there is no tree-level coupling of a neutral scalar or 
vector boson to bd' to be renormalized. Hence we cannot absorb the anti-hermitian parts 
of the flavor-changing self-energy matrices into counterterms and they do contribute to 
our calculation. 

The absorption of the hermitian parts into wave function counterterm matrices is op- 
tional, because the introduction of wave function counterterms only trivially shuffles 
self-energy contributions into vertex counterterms. It is most straightforward then to 
avoid the issue of counterterms altogether by simply calculating the fermion self-energy 
diagrams as they are shown in fig. 0. This calculation is straightforward because the 
internal b quark line is off-shell and therefore it does not contribute to the 1 -particle 
pole of the s quark and needs not be truncated. It is crucial to note that one must start 
with m s 7^ m b in the diagrams with external self-energies in fig. § to properly account 
for the quark propagator l/(m& — m s ), and then take the limit m^m s — > (except in 
the tan /3-enhanced Higgs couplings, of course) at the end. FCNC transitions become 
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meaningless for degenerate quark masses, and one obtains an incorrect result if one 
starts with m s = m b and regulates the propagator pole with an off- shell momentum p 
with p — > 0. In this respect the Higgs exchange diagrams in fig. ^ differ from the situa- 
tion with 7- or Z- exchanges, where both methods give the same correct result. Further 
we note that for m s 7^ m b one must include flavor-changing self-energies in external 
lines with a factor of 1 rather than of 1/2 as in the flavor-conserving case. This is due 
to the fact that flavor-conserving self-energies come from the residue of the one-particle 
pole, while in our approach flavor-changing self-energies are part of the non-truncated 
Green's function. By close inspection of the formulae in JT7[ | one also recovers this 
"factor of 1 rule" from the expressions for the wave function renormalization matrices 
derived in [0~7|]. 

There are two fermion self-energy diagrams that contribute in the 2HDM, one with a 
would-be Goldstone boson G + and one with the physical charged Higgs H + . Their sum 
is ultraviolet-finite. This is different from the SM case, in which the UV divergence of 
the G + diagram cancels with the UV divergence of a SM Higgs vertex diagram involving 
a G + and a top quark in the loop. As in the penguin diagrams involving H + and G + 
in the loop, only the internal top quark contributions to the self-energy are relevant 
here, because the coupling of H + or G + to quarks is proportional to either of the quark 
masses and we neglect m s . The self-energy diagrams add the following term to the 
Wilson coefficients: 



C seif = ^£ tan 2^ ^ Xh+ _ p + ( XH+fXt ) 



sin 2 a cos 2 a 



+ 



M 2 M 2 



2 ™ ~ v ~" ''——Ml 



3.4 2HDM Wilson coefficients and branching ratios 

Adding (jig), Q and we obtain the 2HDM Wilson coefficients in < 



Cs = Cp = r tan p , (16) 

2M^ H r - 1' 

where r = Xg+fxt — Mjj + /m^{m t ). Note that the dependence on the masses of the 
neutral Higgs bosons from the penguin and fermion self-energy diagrams has dropped 
out in their sum without invoking any relation between the mixing angle a and the Higgs 
masses. The result depends on only two of the 2HDM parameters: tan (3 and M H +. 
The two hadronic matrix elements involved are related by the field equation of motion 



(0\hV(x)\B d ,(P Bdl )) = if Bdl PLe 



-iP Bj 
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;0|6 7 fed'(ar)|S <( r(P B J> 



B 4f_ e -iiV x . 



m 6 + m d / 



(17) 



The resulting decay amplitude is 



v^2 47r sin 2 6 1 



(m Bd/ C 5 ) & + ( m Bd ,C P - 



m B. 



(18) 



Here /# is the B^' decay constant, normalized according to = 132 MeV. The cor- 
responding branching ratio is 



B(B d , 



r 2 rv 2 



™?B d ,TB d ,P Bd , 



32vr 2 sin 4 9 
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(19) 



where t b ,, is the B d i lifetime. 

Numerically, the branching fractions are given by 
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•(20) 



It is a well known property of the 2HDM that there exists a limit in which the particles 
A , H°, and H + become very heavy and decouple from processes at the electroweak 
energy scale while h° remains light and its couplings approach those of the SM Higgs 
particle Jl9|]. In the limit of large M H +, Cp and C$ fall as M H + . Thus the deviation of 
the branching fractions from their SM prediction falls as M^ 2 in the large M H + limit, 
and the effects of the enlarged Higgs sector decouple. 

Next we discuss the accuracy of the large tan (3 approximation. Subleading terms in 
tan/5 could be enhanced by powers of m t /mb compared to our result in (|16|), as is 
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the case for the SM contribution. Such terms indeed occur, but they are suppressed 
by two powers of cot (3 compared to the SM terms in (|16|). Hence the formulae above 
are sufficient for all purposes; e.g. if tan/3 = 50 the terms subleading in tan (3 give a 
correction only of 0(2%). If tan (3 is between a few and 15 the 2HDM corrections are 
small and experimentally hard to resolve, so that an error of order cot (3 is tolerable as 
well. 

4 Comparison with other calculations 

4.1 The analyses of He et al. and of Savage 

In the paper of He, Nguyen and Volkas [[3j], the decays B — > £ + £~ , B — > Kl + l~ and 
b — > s£ + £~ are analyzed in both type-I and type-II 2HDMs. In [Tip, the only diagrams 
considered are the box diagram involving two charged Higgs bosons and the A penguin 
involving H + and W + in the loop. Although the calculations of [j^] are performed in 
the 't Hooft-Feynman gauge, the A penguin involving H + and G + in the loop is not 
considered. Similarly, in the paper of Savage [Q], the decay B — > is considered in 

the 2HDM, with and without tree-level flavor-changing neutral Higgs couplings. Only 
the contribution of the A penguin is considered. In both [TJ] and ffi\, several diagrams 
that are important at large tan (3 and required in order to obtain a gauge-independent 
result are neglected. 

4.2 The analysis of Skiba and Kalinowski 

In the paper of Skiba and Kalinowski [§], the decay B s — > r + r~ is analyzed in the 
type-II 2HDM. In [|5|] additional penguin diagrams are considered that are not propor- 
tional to tan 2 (3, but rather contain one or no powers of tan (3. We have neglected these 
contributions in our analysis, because they are not relevant for the interesting case of 
large tan (3. These additional penguin diagrams can be important for small values of 
tan (3 and in regions of the parameter space where some of the Higgs quartic couplings 
are very large resulting in large trilinear H + H~H couplings (H = h°, H°, A ). 
Considering only terms proportional to tan 2 (3, our results for the individual diagrams 
agree with those of with two important exceptions. First, the authors of 0] incor- 
rectly conclude that the box diagram involving H + and W + is negligible and therefore 
neglect it. If we neglect the box diagram, we find that the sum of the remaining con- 
tributions proportional to tan 2 (3 is gauge-dependent. In the 't Hooft-Feynman gauge 
employed in [|5]] the omitted diagram gives the dominant contribution, affecting the nu- 
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merical result substantially. Second, our result for the A penguin diagram involving 
H + and G + in the loop differs from that of [Q] by a sign. Our sign is required for the 
gauge-independence of Cp. 

4.3 The analyses of Huang et al. and Choudhury et al. 

In the paper of Dai, Huang and Huang [pj], the Wilson coefficients in @ are calculated in 
the type-II 2HDM at large tan (3. As in our calculation, only the diagrams proportional 
to tan 2 (3 are considered. However, the authors of [^|] consider only the penguin and 
fermion self-energy diagrams with neutral Higgs boson exchange and neglect the box 
diagram with a W and charged Higgs boson in the loop. Still, after leaving out the box 
diagram, our results for Cs and Cp in the 't Hooft-Feynman gauge do not agree with 
those of [§]. This is partly due to a typographical error in [§], which is corrected in [0,|8p. 
There are two remaining discrepancies. First, our result for the A penguin diagram 
involving H + and W + in the loop differs from that of [§] by a sign. Second, in |Q] a 
contribution from the h° and H° penguin diagrams with two H + bosons in the loop is 
included. This diagram is included in [Bp because it apparently receives a factor of tan (5 
from the trilinear H + H~ H couplings (H = h°, H°). We argue in appendix [A] that the 
trilinear couplings should not be considered tan (5 enhanced. Therefore we conclude 
that the penguin diagram with two H + bosons in the loop should not be included in the 
C(tan 2 (5) calculation because it is of subleading order in tan (3. 

In [0, P the Wilson coefficients in (||) are calculated for supersymmetric models with 
large tan (3. If the diagrams involving supersymmetric particles are neglected, this cal- 
culation reduces to that for the type-II 2HDM with parameters constrained by supersym- 
metric relations. Again, in [0,|8[] only the diagrams with neutral Higgs boson exchange 
are considered, and the box diagram with a charged Higgs boson and W boson is not in- 
cluded. Leaving out the box diagram, our result for Cs and Cp in the 't Hooft-Feynman 
gauge does not agree with the non-SUSY part of that of [0,0] . This discrepancy arises 
because our result for the penguin diagrams involving H + and W + in the loop differs 
from that of [[7|,|8p by a sign. Once SUSY relations are imposed on the Higgs sector, it is 
clear that the penguin diagrams with two H + bosons in the loop are not of order tan 2 (5, 
and the authors of [0,^]] have omitted these diagrams, as we did. 

A final critical remark concerns the treatment of the renormalization group in the paper 
by Choudhury and Gaur [§]. They include an additional renormalization group factor 
to account for the running of the Wilson coefficients. Yet these authors have overlooked 
that the running of the (chiral) scalar quark current in Q s and Q P (see ([7])) is compen- 
sated by the running of the 6-quark mass multiplying the currents as explained in sect. 
This leads to an underestimate of the Wilson coefficients by roughly 23%. 
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In conclusion, the papers in [|5|] and disagree with each other, and our cal- 

culation does not agree with any of them. None of the results in passes the test of 
gauge-independence and, in our opinion, each contains mistakes. 

5 Phenomenology 

As can be seen from the numerical coefficients in (|2"C|), B(B S — > £ + £~) is significantly 
larger than the corresponding branching fraction for B d . This is due primarily to the 
relative sizes of |Vt s | and \V t d\- As a result, even though the production rate of B s is 
three times smaller than that of B d at the Tevatron, the bounds on the leptonic branching 
fractions of B s are much closer to the SM predictions than those of B d [^(J. For this 
reason we concentrate on the decays of B s . Because of the suppression of the branching 
fractions by mf/m 2 B , clearly the decay to r pairs is the largest of the leptonic branch- 
ing fractions in both the SM and the 2HDM. However, this decay is very difficult to 
reconstruct experimentally (due to the two missing neutrinos), and as a result the ex- 
perimental limits on B decays to r pairs are very weak. Therefore in our numerical 
analysis we focus on the decay B s — > for which the experimental bound is the 

closest to the SM prediction. The best experimental bound comes from CDF pH|], where 
one candidate event for B — ► has been reported; this event was consistent with 

the expected background and lay in the overlapping part of the search windows for B d 
and B s . The corresponding 95% confidence level upper bound on the B s — ► /i + yU~ 
branching fraction is [ |20| ] 

B(B S -> < 2.6 x 1(T 6 (expt). (21) 

The SM prediction for the branching fraction is 

B(B S -> /i + /i") = 4.1 x 1(T 9 (SM), (22) 

where we have taken the central values for all inputs in ( |20| ) and ignored the 2HDM 
contributions as well as the errors in the hadronic parameters. 

Because the 2HDM Wilson coefficients in ( |16[ ) depend on only two of the parameters 
of the 2HDM, tan (3 and M H +, the behavior of the result in different parts of parameter 
space is easy to understand. In regions of the parameter space with a large 2HDM 
contribution to B{B d i — > £ + £~) compared to the SM contribution, we may neglect Y{x t ) 
in (PD|). Then the result is particularly simple: the branching fractions are proportional 
to tan 4 /3 log 2 r/(r - l) 2 . 

We can see from ( PD[ ) and the value of Y(x t ) given in (Q) that the interference between 
the SM and 2HDM contributions to the branching fractions is destructive. The effect of 
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Figure 3: B(B S -> in the 2HDM as a function of M H + for tan /3 = 100, 75, 60, 

50 and 25. For comparison we show the current experimental bound [ ]2"0| ] and the SM 
prediction for the branching fraction. The vertical line is the lower bound on M H + in 
the type-H 2HDM from b -»• S7 [gTp. 
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the destructive interference can clearly be seen in fig. [3|. In fig. |3] we plot the predicted 
value of B(B S — ► in the 2HDM as a function of M H +, for various values of tan 8. 

For comparison we show the constraint on M H + from b — > [|ZT|], obtained from the 
current 95% confidence level experimental upper bound of B(b — ► S7) < 4.5 x 10~ 4 
from the CLEO experiment. For very large tan 8 and relatively light H + , the 2HDM 
contribution dominates and the branching fraction is significantly enhanced compared 
to its SM value. As the 2HDM contribution becomes smaller due to increasing M H + 
or decreasing tan (3, the branching fraction drops, eventually falling below the SM 
prediction due to the destructive interference. If we ignore the kinematical factor of 
(1 — Amj/m 2 Bd/ ) in front of C s in ( |T9D (which is a good approximation for B d > — > 
but not for B$ — > t + t~) then we can easily show that the branching fraction in the 
2HDM crosses the SM value when 

m| s tan 2 /3 logr 

8M 2 , — = Y ^ (23) 

and reaches a minimum of half the SM value when 

777.R tan 2 8 logr 1 . , 
Bs , S = -y(x t ). (24) 
8M 2 , r-1 2 v ; v ' 

These correspond to tan 2 (3 log r/(r — 1) = 1790 and 893, respectively. As a numer- 
ical example, taking tan/? = 60 and M H+ = 175 GeV (500 GeV), we find B(B S -> 

= 1-8 x 10~ 8 (2.1 x 10~ 9 ). 
In fig. ^ we plot the regions of M H + and tan (3 parameter space that will be probed as 
the sensitivity to the decay B s — > improves at the Tevatron Run II. The contours 

shown (from left to right) were chosen as follows. The current upper bound on B(B S — ► 



is 2.6 x 10~ 6 from CDF [ ]20| ] with about 100 pb _1 of integrated luminosity. 
This bound excludes a small region of parameter space with very high tan 8 and very 
light H + , shown by the solid line at the far left of fig. [|. Such low H + masses are 
already excluded by the constraint from b — > S7 [|Zl"]]. The rest of the contours in fig. 
^ show the regions that we expect to be probed at the Tevatron Run II and extended 
Run II with various amounts of integrated luminosity, assuming that the background 
for this process remains negligible. In this case the sensitivity to the branching fraction 
should scale with the luminosity. If there is background however, then the sensitivity 
will scale only with the square root of the luminosity. With 2 fb _1 from each of the two 
detectors, the sensitivity should improve by a factor of 40 over the current sensitivity, 
to 6.5 x 10~ 8 , shown by the short dashes in fig. For the values of tan/? that we 
consider, this sensitivity will still only probe values of M H + already excluded by b — > 
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Figure 4: Regions of M H + and tan (3 parameter space probed by measurements of 
B(B S — > Contours are chosen as described in the text. The vertical dashed 

line is the present lower bound on M H + in the type-II 2HDM from b — > 57 [|2T|] . 
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57. We also show two contours for the expected sensitivity with 10 fb _1 and 20 fb _1 of 
integrated luminosity per detector (dotted and dot-dashed lines in fig. 0, respectively). 
These correspond to an extended Run II of the Tevatron. With 10 fb 1 we expect the 
sensitivity to reach a branching fraction of 1.3 x 1CT 8 , allowing one to begin to probe H + 
masses above the current bound from b — > S7 for tan (3 > 54. With 20 fb -1 we expect a 
reach of 6.5 x 1CT 9 , less than a factor of two above the predicted SM branching fraction. 
This would allow one to probe H + masses above the current bound from b — > 57 for 
tan (3 > 47. In particular, for tan (3 = 60, a non-observation of B s — > / u + at this 
sensitivity would rule out H + lighter than 260 GeV. 

Looking farther into the future, the experiments at the CERN LHC expect to observe the 
following numbers of signal (background) events for B s — > fi + \i~ after three years of 
running at low luminosity, assuming the SM branching fraction [J22"|]: ATLAS: 27 (93); 
CMS: 21 (3); and LHCb: 33 (10). Since the suppression of this branching fraction in the 
2HDM is at most a factor of two, the LHC experiments will be able to observe this decay 
for any configuration of the 2HDM at large tan (3. For e.g. tan (3 = 60 and M H + < 285 
GeV, the branching fraction in the 2HDM is enhanced by 30% or more compared to the 
SM. Similarly, for tan (3 = 60 and 375 GeV < M H + < 1 TeV, the branching fraction 
is suppressed by 30% or more compared to the SMP|. In these regions, we expect the 
LHC to be able to distinguish the 2HDM from the SM. In the region of large M H +, 
the dependence on M H + is very weak; hence the LHC measurement will give powerful 
constraints on tan (3 in the large tan (3 region. 

We have made no attempt to simulate the experimental background for this decay in 
order to obtain an accurate estimate of the reach of the Tevatron Run II. Neither have 
we taken into account the theoretical uncertainty. We expect the largest theoretical 
uncertainty to come from uncertainties in the input parameters, primarily the B meson 
decay constants and CKM matrix elements in (|2T)D. These uncertainties will be reduced 
as the B physics experiments progress and lattice calculations improve. Also QCD 
corrections to the 2HDM contribution will arise at NLO and require the calculation of 
two-loop diagrams. In the SM, the NLO corrections increase the decay amplitude by 
about 3%, and therefore increase the branching fraction by about 6%. We expect the 
NLO corrections to the 2HDM contribution to be of the same order, in which case our 
conclusions are not significantly modified. 

In order to evaluate the usefulness of B s — > //+ //~ as a probe of the 2HDM, we must 
compare it to other measurements that constrain the 2HDM in the large tan (3 regime. 
As the statistics of B physics experiments improve, the measurement of b — > S7 will 
improve as well. If B s — ► is to be a useful probe of the 2HDM, it must be sen- 



3 We do not consider charged Higgs masses larger than 1 TeV for naturalness reasons. 
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sitive to a range of parameter space not already explored by b — > S7 at each integrated 
luminosity. Fortunately, B s — ► is complementary to 6 — ► S7 because of the differ- 

ent dependence on tan (3. While the 2HDM contributions to b — > S7 are independent of 
tan /? for tan (3 larger than a few, the 2HDM contributions to B s — > depend sen- 

sitively on tan (3. This makes B s — > an especially sensitive probe of the 2HDM 

in the large tan (3 regime, while b — > 57 will remain more sensitive for moderate and 
small tan (3. Finally, a fit to the Z decay data in the 2HDM [ ]2~3| ] puts weak constraints 
on the H + mass for very large tan/3: M H + > 40 GeV at 95% confidence level for 
tan (3 = 100. The fit gives no constraint for tan (3 < 94. 



6 Conclusions 

In this paper we have analyzed the decays B$ — ► £ + £~ in the type-II 2HDM with large 
tan (3. Although these decays have been studied in a 2HDM before, the previous analy- 
ses omitted the box diagram involving W and H + , which is the dominant contribution 
at large tan (3 in the 't Hooft-Feynman gauge and is needed for gauge independence. We 
showed that when all the contributions are properly included in the large tan (3 limit, the 
resulting expressions for the branching fractions are quite simple and depend only on 
tan (3 and the charged Higgs mass. These 2HDM contributions can enhance or suppress 
the branching fractions by a significant amount compared to their SM values, provid- 
ing tantalizing search possibilities with the potential to probe large parts of the large 
tan (3 parameter space of the 2HDM. We have focused in our numerical analysis on 
B s — > for which the the experimental sensitivity is best. We find that for rea- 

sonable values of tan (3 up to 60 and charged Higgs masses above the lower bound set 
by b — ► S7, B(B S — > can be increased by up to a factor of five above its SM 

expectation or suppressed by up to a factor of two, depending on the charged Higgs 
mass. Although very high statistics will be needed to observe this decay, it promises to 
be an experimentally and theoretically clean probe of new Higgs physics. 
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A Trilinear Higgs couplings 

The trilinear H + H~H couplings (H = h°,H°) for the non-supersymmetric 2HDM 
were first presented in [Q]. These couplings depend strongly on the model consid- 
ered. For the most general CP-conserving 2HDM with natural flavor conservation, the 
H + H~H couplings are given by igQn/M w , where 
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Qh» = — v SM [sin(/3 — a)(2A 3 + A 4 ) — sin/5 cos/3 cos(a + (3)\ 
+2 sin (3 cos (3 (— sin a sin (3\\ + cos a cos /3A 2 )] 

Q H o = — v 2 SM [cos(/3 — a)(2A 3 + A 4 ) — sin /3cos/3sin(a + /3)A 5 

+2 sin j3 cos (3 (cos a sin (3\i + sin a cos /3A 2 )] . (25) 

The H + H~A° coupling is zero. Here vsm = 174 GeV is the SM Higgs vev, and the 
A, are the scalar quartic couplings of the Higgs potential given in To write these 
couplings in terms of Higgs masses and mixing angles, one must make an assumption 
to eliminate one of the independent Aj. In [Q], formulae are presented for the two cases 



A x = A 2 and A 5 = A 6 . The formulae in ([25|) agree with [Q] in these two cases. 
At large tan (3, ( |2"5| ) reduces to 

Q h0 ~ -u| M cosa(2A 3 + A 4 ) [l + 0(cot/3)] 
Q H o ~ -^ M sina(2A 3 + A 4 ) [1 + C(cot /?)]. (26) 
These couplings are not enhanced at large tan (3. 

Considering instead the case Ai = A 2 and writing the trilinear couplings in terms of 
Higgs masses and mixing angles, we find at large tan (3, 

Q h o ~ ~\M 2 H o tan /3 cos 2 a sin a [1 + 0(cot/3)] 
Q H o ~ |M^ tan/3cosasin 2 a [1 + C(cot/5)]. (27) 

Naively, one would conclude that these couplings are enhanced at large tan (3. This is 
incorrect because the angle a depends on tan (3. At large tan (3 we have 

tan 2a = 2 ( 4As + X ^ C ot f3 [1 + 0(cot 2 (3)}. (28) 
A5 — 4(A 2 + A3) 



Thus for generic values of the A i5 sin a ~ cot (3, and the tan (3 enhancement in ( |2"7| ) is 
cancelled. 
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